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ABSTRACT 

This paper reports Very Large Array observations at 325 and 1425 MHz (A90cm and A20cm) 
during and near the periastron passage of HD 80606b on HJD 2454424.86 (2007 November 20). 
We obtain flux density limits (3a) of 1.7 mJy and 48 /jJy at 325 and 1425 MHz, respectively, 
equivalent to planetary luminosity limits of 2.3 x 10 24 erg s _1 and 2.7 x 10 23 erg s _1 . Unfor- 
tunately, these are several orders of magnitude above the nominal Jovian value (at 40 MHz) of 
2 x 10 18 erg s _1 . The motivation for these observations was that the planetary magnetospheric 
emission is driven by a stellar wind-planetary magnetosphere interaction so that the planetary lu- 
minosity would be elevated near periastron. We estimate that, near periastron, HD 80606b might 
be as much as 3000 times more luminous than Jupiter. Recent transit observations of HD 80606b 
provide reasonably stringent constraints on the planetary mass and radius, and, because of the 
planet's highly eccentric orbit, its rotation period is likely to be "pseudo-synchronized" to its 
orbital period, allowing a robust estimate of the former. Consequently, we are able to make 
relatively robust estimates of the emission frequency of the planetary magnetospheric emission 
and find it to be around 60-90 MHz. While this is too low for our reported observations, we 
compare HD 80606b to other high-eccentricity systems and assess the detection possibilities for 
both near-term and more distant future systems. Of the known high eccentricity planets, only 
HD 80606b is likely to be detectable, as the others (HD 20782B and HD 4113) are both lower mass 
and longer rotational periods, which imply weaker magnetic field strengths. We find that both 
the forthcoming "EVLA low band" system, which will operate as low as 65 MHz, and the Low 
Frequency Array may be able to improve upon our planetary luminosity limits for HD 80606b, 
and do so at a more optimum frequency. If the low-frequency component of the Square Kilo- 
metre Array (SKA-lo) and a future lunar radio array are able to approach their thermal noise 
limits, they should be able to detect an HD 80606b-like planet, unless the amount by which the 
planet's luminosity increases is substantially less than the factor of 3000 that we estimate; for the 
SKA-lo, which is to be located in the southern hemisphere, future planetary surveys will have to 
find southern hemisphere equivalents of HD 80606b. 

Subject headings: planetary systems — planets and satellites: magnetic fields — planets and satellites 
(HD 80606b) — radio continuum: planetary systemb 



1. Introduction 

HD 80606 is a G5 star, at a distance of 58.4 pc, 
which is distinguished by being orbited by a planet 
with one of the highest known orbital eccentrici- 
ties. The planet is a giant planet (3.94 ± 0.11 Mj) 
with a 111-day orbital period on an ext r emely 
eccentric orbit (e = 0.93) (jNaef et al.l l200ll: 
Pont et al.|[2009h . 

All of the giant planets in the solar system and 
the Earth generate radio emission as a result of an 
interaction between the solar wind and the plan- 
etary magnctosphcrc. Specifically, the solar wind 
impinging on the planetary magnetosphere gener- 
ates currents within the magnetosphere which are 
then directed into the planet's magnetic polar re- 
gions where an electron cyclotron maser instability 
is produced. Motivated by the radio emission from 
solar system planets, there have been a number of 
searches for magnetospherically generated radio 
emission from extrasolar planets, a t frequencies 



rangi n g from 25 MHz to 14 MHz (lYantis et al 
19771: iwinglee et all !l98i iBastian et al.l 12 000; 



Zarka et al.lboOlULazio et al.ll2004URvabov et al 



2004 1 : iGeorge fc Stevens! l2007t lLazio fe Farrell 
20071 : ISmith et al.l2009tlLecavelier Pes Etangs et al 



20091 ; Lazio et all 2010f) . Sensitivities (la) have 



ranged from roughly 1000 mJy to below 1 mJy, 
with generally better sensitivities obtained at 
the higher frequencies. Typically, these sensi- 
tivities have not been sufficient to challenge the 
predicted levels of extrasolar planetary magneto- 
spheric emission, though some rec ent multi-epoch 
74 MHz observations of t Boo ( Lazio fe Farrelll 



20071 ) do indicate that its planet's emission is not 
consistent with the most optimistic predictions, 
unless the emission is either at a frequency well 
below 74 MHz or is much more strongly beamed 
than that of solar system planets. 

Most recent radio observations of planets have 
focused on planets in orbits with relatively small 
semi- major axes and low eccentricities (but see 



Lazio et al. 2010h . Part of the motivation for this 



focus is the stellar wind powering of the electron 
cyclotron maser instability. Close to the host star, 
the stellar wind should be more intense with a con- 
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comitant increase in the strength of the planetary 
radio emission. A secondary cause is the result of 
a selection effect. Planets close to their host stars 
produce the largest radial velocity signature, and 
the vast majority of known extrasolar planets have 
been found with the radial velocity method. 

A concern with this focus on so-called hot 
Jupiters is that they are likely to be tidally locked 
to their host star, with rotation periods of order a 
few days. To the extent that rapid rotation may 
be necessary for the generation of a strong plane- 
tary magnetic field, these "hot Jupiters" may have 
weak fields that do not generate much electron 
cyclotron maser emission or generate it at a fre- 
quency well below that accessible from the ground. 
Absent a direct radio detection of magnetically 
generated emission, evidence for extrasolar plane- 
tary magnetic fields is ambigu ous. Consiste n t with 
the scenario of weak fields, iFrance et al. ( 201dh 
were unable to find any unambiguous signal of 
aurorally generated UV emission from the "hot 
Jupiter" HD 209458b. From this they derive a 
limit on the intrinsic planetary magnetic moment 
of 0.01 Mj, but perhaps as large as 0.1 TWj, if 
the interaction with the stellar magnetic field is 
important. In turn, this limit implies a mag- 
netic field strength o f only 10% that of Jupiter. 



Batvgin fc Stevenson! (|2010h argue that Ohmic 



dissipation could explain the inflated atmospheres 
of "hot Jupiters," and they are able to explain the 
atmospheres of HD 209458b, HD 189733b, and 
Tres-4b assuming a field strength of 10 G at the 
cloud tops, as compared to the value of 4.2 G 
for Jupiter. Strictly, their estimate (~ 10 G) is 
an upper limit because Ohmic dissipation need 
not be the only effec t operating to i nflate "hot 
Jupiter" atmospheres. lLai et al. I (|2010h show that 
the asymmetric transit properties of WASP-12b 
might be able to be explained by a magnetopause, 
if the planetary magnetic field strength is of order 
1 G, but their model also has a fairly weak depen- 
dence on the magnetic field strength (oc B 1 / 3 ). 
Conversely, a small number of stars with "hot 
Jupiters" s how evidence for sta r-planet interac- 
tions (e.g., IShkolnik et al 



2008) 



which appear 
to have a magnetic component. In contrast to 
"hot Jupiters," with its large eccentricity orbit, 
HD 80606b may be immune to some of the ef- 
fects that might act to damp the magnetic field 
strength of "hot Jupiters." 



2 



The planet orbiting HD 80606 passed through 
periastron on HJD 2454424.86 (2007 November 20 
at 08:05 UT). This paper reports observations 
at 325 and 1425 MHz (A92 cm and A21 cm, re- 
spectively) of HD 80606b close to its periastron. 
In Section [2] we describe the observations, analy- 
ses, and present the images obtained; in Section [3] 
we discuss our upper limits in the context of re- 
cent multi- wavelength observations of HD 80606b, 
and in Section U we make predictions about future 
radio observations. 

2. Observations 

HD 80606 was observed^ at 325 and 1425 MHz 
with the Very Large Array (VLA), in its B config- 
uration, near the periastron passage of the planet. 
Table [1] summarizes relevant details of the ob- 
servations, with specific information about the 
actual times of observations and the planetary 
phase. In addition to HD 80606, at both fre- 
quencies, 3C 286 was observed to serve as a flux 
density and bandpass calibrator, and the VLA 
calibrator J0929+502 was observed as a visibility 
phase calibrator . Although the VL A could oper- 
ate at 74 MHz (jKassim et al.ll2007l ). the 74 MHz 
dipoles were not mounted at the time of these ob- 
servations. 

Standard data reduction techniques were used 
within the AIPSH software package to calibrate 
the spectral bandpass, the visibility amplitudes, 
and visibility phases. The data were acquired in 
a spectral-line mode both to enable the identifica- 
tion and excision of radio frequency interference 
(RFI) and for the purposes of wide-field imaging. 
These issues of RFI excision and wide-field imag- 
ing are particularly relevant for the 325 MHz ob- 
servations. 

Following the data calibration, several itera- 
tions of imaging and self-calibration were used to 
produce images of the field around HD 80606; 
at both frequencies, polyhedral imaging was 
used to image the prim ary beam field of view 
(|Cornwell & PerlevMl992l ). Figures[T]and[2]present 
the inner portions of the resulting images, centered 
on HD 80606. The resulting image noise levels are 
16 //Jy beam -1 at 1425 MHz and 580 /iJy beam" 1 




RIGHT ASCENSION (J2000) 



Fig. 1.— Field around HD 80606 at 1425 MHz 
on HJD 2454425, which includes the time of pe- 
riastron (Table [T]). The cross marks the loca- 
tion of the star. The beam for this image is 
5'.'49 x 4'.'35, the noise level in the image is 
16 //Jy beam -1 , and the gray scale is logarithmic 
between — 80^Jybeam _1 and 1.6 mJy beam -1 . 
The source to the southeast of the star is 
FIRST J092239.6+503529. 



at 325 MHz, which in both cases are within ap- 
proximately 50% of the expected thermal noise 
levels. 

3. Discussion 

At neither frequency do we detect statistically 
significant emission from the direction of the star. 
Table [2] summarizes the resulting flux density lim- 
its on this system near periastron passage and 
the resulting luminosity limits. For reference, the 
nominal luminosity of Jupiter is 2 x 10 18 erg s" 1 
at a frequency near 40 MHz. 

The age of HD 80606 is uncertain, with es- 
timates ranging f rom as little as 0.3 Gyr to in 
exces s of 14 Gvr dWright et al.ll2004 ISaffe et al ' 



2005; IValenti fc Fischeril2005l iTakeda, et all [2007; 



1 program AL715. 

2 version 31DEC08. 



iHolmberg et all l2009t ) ." but with many estimates 
favoring an age of approximately 7 Gyr. For the 
purposes of this discussion, we take the age of 
HD 80606 to be comparable to that of the Sun, 
with a comparable stellar wind strength. Con- 
sequently, we assume that a Jovian planet, with 
an orbital semi-major axis comparable to that of 
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Table 1 
Observation Log 



Frequency 


Wavelength 


Bandwidth 


Epoch 


Orbital Phase 


325 MHz 


90 cm 


12.5 MHz 


HJD 2454423.89-2454424.22 


0.69-0.73 








(2007 November 19 09:20-17:17 UT) 




1425 MHz 


20 cm 


50 MHz 


HJD 2454424.74-2454425.07 


0.94-0.12 








(2007 November 20 05:46-13:44 UT) 





Table 2 

Luminosity Limits for HD 80606b 



Wavelength 


Flux Density Limit 


Luminosity Limit 




(3a) 


(erg s" 1 ) 


1425 MHz 


48 fiJy 


2.7 x 10 23 


325 MHz 


1.7 mjy 


2.3 x 10 24 



Note. — In converting from flux density to luminosity, 
an emission bandwidth of order the observing frequency is 
assumed. 



Jupiter (« 5 AU), would have a radio luminosity 
similar to that of Jupiter. 

It is clear that our luminosity limits remain sev- 
eral orders of magnitude (10 5 -10 6 ) above the lu- 
minosity of Jupiter. The original motivation for 
these observations was the expectation that the 
high eccentricity of the planet's orbit could pro- 
duce a dramatic variation in the strength of the 
planet's radio emission. Because a planet's mag- 
netospheric emission is controlled by the input 
stellar wind power, the luminosity of the planet 
is a strong function of its distance from the star. 
Farrell et al.l ( 19991 Appendix A) discuss the scal- 
ing of the luminosity with distance and find that 
an approximate scaling la w is d x , with x « —1.6. 
Griessmeier et al. I <l2007bllah consider more com- 
plex stellar wind models, but a strong dependence 
on the planet-star distance remains. 

Over the course of the planet's orbit, from apas- 
tron to periastron, we therefore expect that the 
luminosity of HD 80606b could vary by a factor 
of 200. Even at apastron (0.87 AU), HD 80606b 
is still closer to HD 80606 than Jupiter is to the 



Sun. If we consider the further boost that could 
be obtained by moving a Jovian-like planet from 
a distance of 5 AU to 0.87 AU, the boost would 
be a factor of 16. Thus, one might expect that, 
near periastron, HD 80606b could be as much as a 
factor of 3000 times more luminous than Jupiter. 
However, even an increase by a factor of order 3000 
over Jupiter's luminosity would be insufficient to 
increase HD 80606b to a level comparable to our 
detection thresholds. 

A simple explanation for our non-detections is 
therefore that, even with the added boost due to 
the small planet-star distance during periastron, 
the planet's luminosity remains too low to be de- 
tected with current equipment. There are several 
other factors, however, that could also explain our 
non-detection. 

One possibility is that the radio emission is 
not bea med in our directi on. In the case of 
Jupiter, IZarka et al.l ( 2004 ) determined that it 



emits into a solid angle of approximately 1.6 sr, 
directed roughly along Jupiter's rotational axis. 
HD 80606b is one of the rapidly growing class of 
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RIGHT ASCENSION (J2000) 

Fig. 2.— Field around HD 80606 at 325 MHz 
on HJD 2454424, just before the planet's pe- 
riastron passage (Table [1]). The cross marks 
the location of the star. The beam for this 
image is 22'.' 1 x 20'.' 1 and is shown in the 
lower left corner; the noise level in the image is 
580 jitJy beam -1 ; and the gray scale is logarithmic 
between — 2.3 mJy beam -1 and 58 mJy beam -1 . 
The area shown is approximately the same size as 
that shown in Figure [TJ 

transiting extrasolar planets, and recent transit 
observations demonstrate that there is a signifi- 
cant spin-orbit misalignment (~ 50°) in this sys- 



tem (iFossev et al.l200alGarcia-Melendo k McCullougb 
20091 ilMoutou et al.ll2009l : IPont et al.ll2009l : IWinn et al.l 



2009). For a transiting planet whose rotational 
axis is aligned with the stellar rotation axis and 
no spin-orbit misalignment, the probability of in- 
tercepting the emission beam would be relatively 
small. However, the existing spin-orbit misalign- 
ment in this system increases the probability that 
emission beam intersects the line of sight. Of 
course, it is possible that the emission cone of 
HD 80606b is much smaller than that of Jupiter, 
such that the probability of intercept remains 
quite low. 

A second possibility is that the expected up- 
per frequency limit of emission is (well) below 
our observation frequencies. The intense cyclotron 
maser emission of Jupiter cuts off dramatically at 
frequencies above approximately 40 MHz, as the 
local plasma frequency in the polar emitting re- 
gion exceeds the electron cyclotron frequency and 



the radio emission cannot escape. For most ex- 
trasolar planets, one can estimate the upper fre- 
quency limit by scaling their magnetic moments 
from that of Jupiter. Previous estimates for the 
magnetic moment, and upper frequency limit, for 
HD 80606b have varied wildly. lLazio et al.l (|2004l ) 
estimated an up per frequency lim i t of a round 



(|2007bh esti- 



180 MHz, while iGriessmeier et al 
mated 0.8 MHz. 

The combination of the eccentric orbit and the 
planetary transits, however, now provide sufficient 
information to estimate, with relatively low uncer- 
tainty, the upper emission freq uency. This upper 
emission frequenc y is given by (|Farrell et al.l l"l999: 



Lazio et al . 2004) 



24 MHz 



M \ 



5/3 



R 



(1) 



for a given planetary mass M, rotation and ra- 
dius R. In practice, this value may underestimate 
the actual upper emission frequency, as the case 
for Jupiter is v c « 40 MHz whereas Equation ([1]) 
suggests 24 MHz. This expression also assumes 
that the planetary magnetic dipole is not offset 
significantly (as it is for the planets Uranus and 
Neptune), which can produce significant asymmet- 
ric emission. 

Of the various determinations of the mass 



and radius of the planet (Fossev et al. 2009; 


Garcia-Melendo & McCullouffhl2009: 


Moutou et al. 


20091 Pont et al. 20091: IWinn et al. 


2009b. we 



adopt those that iPont et al.l (J2009J) obtained from 
a Bayesian analysis of all of the data avail- 
able to them. They found a planetary mass of 
3.94 ±0.11 Afj and a radius of 0.98 ± 0.03 Rj, 
implying a density about four times higher than 
that of Jupiter; their models also suggest that the 
core mass of HD 80606b could be several times 
that of Jupiter. 



Hutl (|198lh discusses the tidal effects in eccen- 
tric systems. The large eccentricity for the orbit 
of HD 80606b should have driven it into a state of 
pseudo-synchronization between its rotation and 
orbita l periods. Using expressions derived by [Hut 
(|198ll ). we estimate a (pseudo-synchronized) rota- 
tional period of 39.9 hr. 

Evaluating Equation ([I]), we find that the upper 
emission frequency of approximately 55 MHz, and 
potentially as high as 90 MHz, if one accounts for 
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the 50% difference between the actual value for 
Jupiter and that predicted by Equation (TTJ). It 
is clear that this upper emission frequency is well 
below that of our observations. 

This estimate for the em i ssion frequency dif- 
fers from that of Lazio et al. | (|2004 by only a fac- 
tor of 2-3, but it differs from that obtained by 
iGriessmeier etHl (l2007bh by a factor of order 10 . 
We believe that the differences can be attributed 



to ass umptions of the rotation period. lLazio et al 



([20041 ) assumed a rotation period equal to that of 
Jupiter, which is within a factor of 4 of the rotation 



period we adopt here. However. IGriessmeier et al 



(|2007bt ) identified the planet as "tidally locked," 
which would mean that it would have an extremely 
slow rotation rate, and therefore a low magnetic 
moment and low emission freq uency. One of the 
key fac tors that may have led IGriessmeier et al 



(|2007bf ) to identify the planet as tidally locked is 
that they appear to have used expressions for es- 
timating the synchronization time scale which as- 
sume a circular orbit (their Appendix B). 

4. Future Prospects 

We consider the possibilities of future observa- 
tions from two perspectives. First, are there plan- 
ets other than HD 80606b that should be targeted? 
Second, what are the prospects for future instru- 
ments to detect HD 80606b? 

While, HD 80606b is not the only known 
planet with a high-eccentricity orbit, it appears 
the most promising for future radio observations. 
HD 20782b is a 1.9 Mj (minimum mass) planet 
in an orbit with an eccent ricity of 0.97 and a n 
orbital period of 592 days ((O'Toole et al.ll2009h . 
and HD 4113b is a 1.6 Mj (minimum mass) planet 
in an orbit with an eccentri city of 0.903 and a n 
orbital period of 527 days (jTamuz et al . 2008). 
The pseudo-synchronization rotational period is 
coupled to the orbital period, so these longer or- 
bital periods imply substantially larger rotational 
periods, as compared to HD 80606b. We esti- 
mate rotational periods of 63.8 hr and 14 days 
for HD 20782b and HD 4113b, respectively. Com- 
bined with their lower masses, Equation (fTJ) sug- 
gests that their emission frequencies are likely 
to be below 20 MHz. Unless future surveys re- 
veal planets with masses and orbital eccentrici- 
ties comparable to those of HD 80606b (but see 



O 'Toole et aDl2009h . it appears likely to remain 
the most promising target for radio observations 
of (high-eccentricity) planets. 

With respect to observations of HD 80606b, our 
analysis suggests that future observations must be 
both more sensitive and at a lower frequency than 
those that we report here. There are several ob- 
serving systems that are either under construction 
or in the design phase that we now consider. 

The Expanded Very Large Array (EVLA) 
project exposed ce rtain shortcomings o f the VLA's 
74 MHz system (jKassim et all l2007h . primarily 
its generation of harmonics at higher frequencies 
that fall within the EVLA's operational spectral 
window. Consequently, the NRL has funded the 
National Radio Astronomy Observatory to de- 
velop a new suite of low-frequency receivers. With 
the previous 74 MHz VLA system, automated 
data processing of approximately a 1 hr obser- 
vation could obta in rms intensity le vels of order 
100 mJy beam" 1 (|Cohen et al.ll2007l) . while more 
focused efforts on longer observatio ns could ob- 
tain noise levels of 25 mJy beam -1 (jCohen et al 



2004). At a minimum, the new "EVLA- low band" 



system will have a substantially larger bandwidth. 
There has also been continuing algorithmic im- 
provements in the areas of RFI identifica tion and 
excis i on and ionospheric ca l ibration (e.g., Athreyal 
20091 llntema et all l2009t lOffringa et all l2010h . 



As an ansatz, we assume that a combination of 
algorithmic improvements and the larger band- 
width could produce an approximate factor of 10 
improvement in sensitivity, i.e., 3 mJy beam -1 . 
The resulting (3a) luminosity limit would then 
be 10 23 erg s _1 (for an assumed 10 MHz intrin- 
sic emission bandwidth of the electron cyclotron 
maser emission). This value is comparable to our 
more stringent limit (at 1425 MHz), but it would 
be at a frequency at which emission might more 
likely be detected. 

The Low Frequency Array (LOFAFi[f| ) is cur- 
rently under construction in the Netherlands and 
other European countries. When complete, it 
will consist of many dipole phased-array stations, 
operating in a "high band" (HBA) and a "low 
band" (LBA). LOFAR-LBA band will cover the 
frequency range 30-80 MHz, with a peak sensi- 
tivity around 60 MHz, well matched to the ex- 



3 http: / / www. astron. nl / radio-observatory / astronomers/lofar-astronomers 
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pected emission frequency of HD 80606b. Cur- 
rent estimates for the sensitivity of LOFAR-LBA 
are that it will produce a similar or slightly bet- 
ter thermal noise level (about 7 mJy beam -1 in a 
1 hr integration) as the EVLA-low band system. 
We emphasize that the limiting factor in imag- 
ing sensitivity might very well not be the ther- 
mal noise level, but other factors such as iono- 
spheric calibration. Nonetheless, it appears that 
the LOFAR-LBA band system would produce sim- 
ilar constraints. 

Following the EVLA-low band, LOFAR, and 
other pathfinder systems are likely to be the low- 
frequency component of the Square Kilometre Ar- 
ray (SKA-lo) and ultimately a low frequency ar- 
ray on the Moon (the Lunar Radio Array). The 
current design goal for SKA-lo has a lower op- 
erational frequency of 70 MHz with enough col- 
lecting area that the notional thermal noise level 
would be of order 10 /uJy beam -1 in a 1 hr in- 
tegration. Unfortunately, candidate sites for the 
SKA-lo are in the southern hemisphere (Australia 
and South Africa, terrestrial latitude w —30°), 
from which HD 80606b cannot be observed. If 
future planetary surveys detect other HD 80606b- 
like planets and if the thermal noise level can be 
approached with the SKA-lo, the implied luminos- 
ity limit would then be 3 x 10 20 erg s _1 , well within 
the range we estimate might be plausible based 
on the larger stellar wind loading of HD 80606b's 
magnctosphere. 

We thank G. Laughlin for the initial inspira- 
tion for these observations and for pointing out 
that HD 80606b should be in a state of pseudo- 
synchronization; J. Schneider for the Extrasolar 
Planets Encyclopedia; and P. Perley and W. Koski 
for helpful discussions. T.J.W.L. thanks USNO, 
where much of this paper was written, for its hos- 
pitality. The National Radio Astronomy Obser- 
vatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by 
Associated Universities, Inc. The LUNAR consor- 
tium is funded by the NASA Lunar Science Insti- 
tute (via Cooperative Agreement NNA09DB30A) 
to investigate concepts for astrophysical observa- 
tories on the Moon. Basic research at NRL is sup- 
ported by 6.1 Base funding. 
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